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I.  INTRODUCTION 


The  ability  of  a  monopulse  processor  to  determine  the  angular  direction  of 
an  incident  signal  is  limited  not  only  by  the  inherent  frcnt-end  receiver  noise 
but  also  by  the  effect  of  interfering  signals.  If  there  are  other  time  coincident 
signals  present  that  interfere  with  the  signal  whose  direction  is  to  be  estimated 
then  McAulay  [1]  has  shown  that  bias  effects  occur  that  can  seriously  degrade  the 
quality  of  the  estimate.  It  becomes  important  therefore,  to  know  when  such  interference 
is  present  so  that  low  confidence  can  be  assigned  to  the  associated  azimuth  estimate. 

Sherman  [2]  has  observed  that  when  interfering  signals  are  present  the  out¬ 
puts  of  the  monopulse  sum  and  difference  beams  become  incoherent.  He  proposes 
to  use  the  quadrature  information  to  resolve  the  target  and  interference  signals. 

Our  approach  is  to  develop  an  interference  flag  that  indicates  when  more  than 
one  signal  is  present  In  ue  receiver  channels.  Depending  on  the  application, 
the  flag  would  be  used  to  assign  a  low  confidence  to  tne  associated  angle  esti¬ 
mates  or  to  delete  the  angle  estimate  altogether.  In  Section  II,  we  formulate 
the  test  for  interference  as  a  hypothesis  test.  Using  the  Generalized  Likelihood 
Ratio  test  [3]  and  following  the  analysis  of  Hof s tetter  and  DeLong  [4],  we 
obtain  the  optimum  interference  detection  statistic. 

The  performance  of  the  detector  is  analyzed  in  Section  III  where  it  is  first 
shown  that  the  interference  statistic  has  the  Rician  distribution  Exact  evalu¬ 
ation  of  the  false  alarm  and  detection  probabilities  becomes  intractable  and 
use  is  made  of  the  Gaussian  approximation  to  the  Rician  variate.  Numerical 


1 


tpj V^TS^J  ^^r^CTfi-riP}- ■,  ,  Tf*  xap-  -«!**  J-tU 

;»•;,»••;;  . 


results  for  some  typical  cases  of  interest  are  given  and  it  is  shown  that  the 
results  depend  strongly  on  the  relative  phase  between  the  target  and  interfering 
signals,  but  that  good  overall  performance  can  be  obtained.  Conditions  under 
which  the  Gaussian  approximation  is  valid  are  given  and  are  shown  to  hold  for 
the  cases  studied. 


2 


II.  PROBLEM  FORMULATION  AND  SOLUTION 

We  shall  restrict  our  attention  to  amplitude-comparison  monopulse  processing 
that  is  performed  on  a  sampled-data  basis.  Assuming  mixer  preamplifiers  at  the 
output  of  each  antenna  beam  channel,  the  received  signal  samples  in  the  absence 
of  interference  are  modeled  by 

y.j  ~  A^  e  +  n.  1  -  l,2,...,m  (1) 

where  y^  refer  to  the  complex  output  of  ith  antenna  beam  channel;  As’^S’6S’aS  are 
the  amplitude,  phase,  azimuth,  and  elevation  of  the  target  signal;  G . ( - )  is  the 

j.  h 

antenna  patterns  of  the  i  antenna  beam  which  may  be  complex  in  general; 
represents  zero  mean  Gaussian  noise  samples  due  to  the  mixer  preamplifiers  whose 

p 

real  and  imaginary  parts  have  variance  0  •  As  shown  by  Hofstetter  and  DeLong  [4], 
this  model  arises  when  the  received  signal  is  preprocessed  by  a  matched  filter. 

It  can  also  be  used  to  describe  the  case  In  which  a  simple  on-off  pulse  is  trans¬ 
mitted  and  preprocessed  by  a  filter  whose  bandwidth  is  at  least  equal  to  the 
reciprocal  of  the  rise  time. 

If  interference  is  present  that  overlaps  the  target  signal  return,  the 
received  signal  samples  can  be  written  as 

j<Pc  M 

y1  =  A^  e  Aj  e  G,j(9j>  ^  a  l»2,...»m 

(2) 
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where  now  Aj,  6j,  aj  represent  the  amplitude,  phase,  azimuth,  and  elevation 
angles  of  the  interference  signal.  In  addition  to  describing  the  effects  of 
multipath  in  an  L-band  radar  the  model  also  arises  in  the  design  of  the  Discrete 
Address  Beacon  System  (DABS)  that  is  to  be  used  to  perform  the  surveillance 
function  in  the  next  generation  Air  Traffic  Control  (ATC)  system  [5].  During  the 
transition  from  the  present  Air  Traffic  Control  Radar  Beacon  System  (ATCRBS)  to 
a  completely  DABS  operation,  ATCRBS  will  represent  a  source  of  Interference  to 
DABS.  Since  direction  finding  is  performed  at  L-band  using  simple  on-off  pulses, 
the  received  signal  samples  will  be  described  by  (2).  We  note  that  in  this  case, 
however,  that  the  downlink  carrier  frequency  is  known  only  to  within  +  3  MHz. 
Therefore,  the  preprossing  filter  cannot  be  exactly  matched  to  the  downlink 
signal  and  a  slight  reduction  in  signal-to-noise  ratio  (SNR)  must  be  tolerated. 

Using  (1)  and  (2)  we  can  formulate  the  following  hypothesis  test  for  the 
detection  of  interference: 

Hgi  interference  absent 

y.  =  A$  e  Gi(6s,as)  +  n.,  i=l,...m  (3) 


H.j :  interference  present 


yi  =  AS  6  Gi(0s,as^  +  AI  e  Gi\ei*ai)  +  ni  •  (4) 


The  solution  to  this  hypothesis  testing  problem  can  be  obtained  from  application 
of  the  Generalized  Likelihood  Ratio  criterion  [3].  First,  we  form  the  ratio 
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where  we  have  used  the  notation  ^  =  (y-j  »y2»*-*  »yrn)  to  denote  the  complex  data 
vector  and  £_  =  {A,</>,6,a)  to  denote  the  unknown  parameter  vector,  and  where  p(yjHk) 
denotes  the  probability  density  function  of  ^  under  either  hypothesis  HQ  or  H^. 

Once  A(jr)  is  computed,  interference  is  declared  present  if  and  only  if  A(^)  >  X\ 
Since  the  noise  components  are  Gaussian  random  variables,  then  under  the 
null  hypothesis,  we  can  write 


max  p(;/J  Hq,£s) 
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(6) 

Hofstetter  and  DeLong  [4]  have  shown  that  if  both  es  and  are  unknown  then  for 
reasonable  antenna  patterns,  G.|(es»a$)>  the  maximization  in  (6)  can  be  solved  if 
there  are  3  linearly  independent  antenna  beams.  If  only  es  is  unknown  or  if  G^e^) 
is  a  fan  beam  in  elevation,  then  they  shown  that  (6)  can  be  solved  if  there  are 
2  linearly  independent  antenna  beams.  With  these  restrictions  on  m,  namely 
m=2  or  3,  depending  on  whether  e<.  and/or  a$  are  unknown,  Hofstetter  and  DeLong 

AAA  A 

show  that  (6)  is  maximized  at  the  parameter  values  A,<p,e  and/or  a  where 
in 

Y  6,(8, a)  Re(y1e-J<p) 
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and  hence  A, <p,6, a  denote  the  maximum  likelihood  estimates  of  A,<p,6,a. 
It  is  also  shown  that  the  maximum  value  of  the  density  function  is 


max 
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Under  the  alternate  hypothesis  the  density  function  is 


P(ztH!  =(2tt  o2)’m/2  exp 
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(11) 


When  both  es  and  <*$  are  unknown,  we  require  m  =  3  and  note  that  there  are  6 
measurements  and  6  unknown  parameters.  If  only  6$  is  unknown,  we  require  m  *  2 
and  then  there  are  4  measurements  and  4  unknown  parameters.  In  either  case, 
the  maximum  value  of  the  density  is  (2tt  o2)"^2  and  it  is  achieved  by  picking  the 
parameter  estimates  to  solve  the  equations 
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Using  these  facts,  the  Generalized  Likelihood  Ratio,  (5),  becomes 
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hence,  we  declare  interference  present  (hypothesis  H^)  if  and  only  if 


I  N2-lXy?l>*  • 
1=1  1=1 


In  the  application  to  the  ATC  problem  only  the  aircraft  azimuth  can  be 


estimated  since  elevation  fan  beams  are  used  to  provide  surveillance  for  high 


and  low  altitude  aircraft.  In  this  case,  two  antenna  beams  having  azimuth 


directionality  would  be  used  that  would  take  the  form  of  a  sum  (even)  and  dif¬ 
ference  (odd)  monopulse  configuration.  Letting  y^ ,  y2  denote  the  outputs  of  the 
sum  and  difference  beams  respectively,  then  interference  will  be  declared 


present  if  and  only  if 


1^1 1  +  i  ^2 1  1^1  +  ^2 I  >  x 


Following  [4],  we  recognize  that  an  equivalent  test  to  (15)  is  to  declare 


Interference  absent  if,  and  only  if, 


-  x  <  jyi  +  j  y2j  -  |y-|  “  J  x 


It  Is  interesting  to  note  that  in  the  test  for  signal  plus  noise  versus  noise 
alone,  which  was  the  problem  that  was  considered  in  [4],  the  two  quantities  in 
(16)  are  added  rather  than  subtracted.  Furthermore,  the  maximum  likelihood 
azimuth  estimate  is  related  to  the  phase  difference  between  the  signals  y^  +  j  y2. 
Since  these  signals  are  readily  formed  at  RF,  it  is  possible  to  obtain  the 
detection  statistic,  interference  statistic,  and  azimuth  estimate  from  the  same 
RF  hardware  configuration.  In  the  remainder  of  this  paper  we  shall  limit  our 
attention  only  to  the  case  of  a  two  beam  monopulse  radar,  hence  we  assume  that 
the  antenna  patterns  depend  only  on  one  angular  variable.  For  convenience  we 
shall  deal  with  the  estimate  of  azimuth,  but  it  is  obvious  that  the  results 
apply  directly  to  a  monopulse  system  that  tracks  in  elevation. 


III.  DERIVATION  OF  THE  EXACT  PDF 

Since  our  basic  motivation  for  undertaking  this  problem  lies  in  t? ■:;  ATC 
context,  we  can  assume  that  m  =  2  and  focus  our  attention  on  (16)  Aid  attempt  to 
analyze  the  performance  of  the  detector  by  computing  the  false  alarm  and  detec¬ 
tion  probabilities.  Therefore,  we  shall  all  attempt  to  compute  the  probability 
density  function  (pdf)  of  the  detection  statistic 

t  *  |y]  +  j  y2|  -  |y1  -  j  y2|  •  (17) 

We  begin  the  analysis  by  defining  new  random  variables 


y±  =  Vl  ±  J  y2 


j(£V  J^T 

=  As[G1(es)  ±  j  G2(es)]  e  0  +  Aj[G1  (ez)  ±  j  G^Qj)  e  ]  +  nj  i  i  n2 


^  e^S  G1  (e$)  t-  Aj  e^1  G]  (Oj) 

±  J 

tyc  j(/?T 

AS  e  0  G2(es)  +  Aj  e  G2(ez) 

• 

(; 

where  we  have  defined  r.ew  noise  variables 


2 

which  we  note  are  zero  mean.  Independent  complex  Gaussian  with  variance  4 o. 


Therefore,  y+  are  also  independent  complex  Gaussian  random  variables  with 
variance  4o^  and  means 


- 


J(Pc 


A$  e  G^Sj)  +  Aj  e  G]  (6j) 


±  J 


AS  e^S  G2^eS^  +  AI  ^  W 


(20) 


Therefore  the  random  variables 


(21) 


are  independent  Rician  random  variables  of  order  2,  [6],  and  their  probability 
density  functions  (pdf)  are  therefore 


P{  (u) 
x+ 


(22) 


where  I  (•)  is  the  modified  Bessel  Function  of  the  first  kind  of  order  zero, 
0 

and  where 


i\  -  IwJ2  ■  (23) 

Using  the  definitions  in  (17),  (18),  and  (21)  we  see  that  the  detection  statistic 
is  given  by 
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hence  its  pdf  is 


P£(u)  =  j  p^(u  -  v)  (v)  dv  . 


(25) 


Then  using  (22)  in  (25)  it  can  be  shown  that  the  exact  form  for  this  pdf  is 


v/^a 


pt(u) 


2 

2  a 


f(x,u//20)  dx 


(26) 


u/f2o 


where 


f(x,y)  =  (x2  -  y2)  exp 


IQ[d+(x  +  y)]  IQCd_ (x  -  y)] 

(27) 


Since  it  is  difficult,  if  not  impossible,  to  evaluate  (26)  analytically  or 

numerically,  we  have  found  it  more  appropriate  to  approximate  the  Rician  variate 

by  the  Gaussian  density.  In  Fig.  1  we  have  shown  curves  obtained  from  [7] 

2  2 

that  show  the  Rician  pdf  for  several  values  of  the  parameter  dy2 a  .  From  the 
figure  it  appears  that  for  values  of  d+/ J 2 a  greater  than  3  the  Gaussian  approxi 
mation  is  quite  good.  In  the  Appendix  we  perform  a  detailed  error  analysis 
that  shows  that  except  for  2M  of  the  area  under  the  lower  tai . ,  the  Rician 
pdf  is  well  approximated  by  the  Gaussian  density  provided 


(23) 


Therefore  we  can  write  (22)  as 


\(u) :  -i 


4ir  a 


exp 


(u  -  d+y 


(29) 


from  (24)  the  detection  statistic  is  i  =  Z+  -  so  that  z  is  also  Gaussian 


bu t  with  mean  d+  -  d_  and  variance  4  a  .  Therefore 


p*(u) 


\/8ti  a 


exp 


[u  -  (d+  -  d.)f  | 
1 — 


8  o 


(30) 


where  d^re  given  by  (23)  and  (20). 

In  the  next  section  we  shall  compute  the  false  alarm  and  detection  prob¬ 
abilities  and  hence  develop  the  criteria  needed  to  evaluate  the  performance  of 
the  detector. 
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IV.  FALSE  ALARM  AND  DETECTION  PROBABILITIES 


Since  the  detection  statistic  can  be  well  approximated  by  a  Gaussian  random 
variable,  it  is  a  straightforward  problem  to  calculate  the  false  alarm  and  de¬ 
tection  probabilities. 

a)  False  Alarm  Probability 

In  this  case,  the  interfering  signal  is  absent  hence  Aj  =  0. 
Furthermore,  the  target  of  interest  lies  within  the  mainbeam  of  the  antenna  and 
it  can  therefore  be  assumed  that  G . ( ©s )  is  real.  Using  these  facts  in  (20), 

(23)  becomes 


G^(es)  +  G|(es) 


(31) 


Therefore  the  mean  value  of  the  detection  statistic  is 


l  =  d+  -  d_  =  0  .  (32) 

A  false  alarm  is  made  whenever  |i,|  >  x,  hence  the  false  alarm  probability  is 


=  2  erf(j^)  (33) 
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where 


erf(x) 


.-l/V 

V  2z_L 


2 

-r/2 


dt 


(34) 


This  shows  thac  -he  detection  threshold  can  be  set  once  the  level  of  the 
backgrourJ  noise  is  known. 

To  verify  that  the  Gaussian  approximation  is  indeed  valid  in  this  case,  we 
see  from  (28)  that  it  is  sufficient  to  have 


2  0 
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10  . 


(35) 


2  2  2 

The  detection  sigral-to-noise  ratio  is  A<.  G^(es)/2cr  and  since  this  quantity  is 
at  least  20  dB  in-the  ATC  context  we  see  that  (35)  will  easily  be  satisfied, 
b)  Detection  Probability 

In  this  case,  we  detect  interference  when  |£|  >  x,  hence  if  the 
Gaussian  approximation  is  valid 


v  ■ 
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exp 


-X  \/8ir  a 


[u  -  (d+  -  dj] 
T~~ 


8  c 


du 


=  erf 


A  -  (d+  -  dj 


To" 


+  erf 


A  +  (d+  -  dj 


2o 


(36) 


If  we  desire  a  false  alarm  probability  PpA,  then  we  solve 


PFA  =  2  erf(ApA) 


(37) 
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for  Then  setting  the  threshold  according  to 

A  =  2o  ApA  (38) 

we  see  that  the  detection  probability  is 

PD  ’  erf(JFA  -  -V-1)  +  erff FA  +  -V"1)  '  (39) 

We  have  yet  to  verify  the  validity  of  the  Gaussian  approximation  in  this 
case,  but  this  requires  the  evaluation  of  (23).  We  shall  consider  this  point 
in  detail  in  the  next  section.  At  that  time,  we  shall  consider  specific  inter¬ 
ference  cases  of  interest  and  carry  out  the  evaluation  of  the  receiver  performance 


in  detail . 


V.  RECEIVER  PERFORMANCE 

In  many  cases  of  practical  interest  the  antenna  patterns  will  be  complex  at 
least  in  a  region  beyond  the  near*in  sidelobes.  Since  the  interference  can 
originate  from  any  azimuth,  we  make  the  complex  dependence  clear  by  writing 


Gi (0j)  =  A1  (_6j )  exp[j  ^ ( ex ) ] 


(40) 


Substituting  (40)  into  (20)  and  using  the  fact  that  G..(es)  is  real  we  can  shown 
after  tedious  but  straightforward  manipulations  that 


~ =  j G-j  ( Qg )  +  ^(^SM  +  P  ±  2  A^  A2  sin(^  -  ^)  +  Ag 
AS 


+  2r 


1  /2 

(0$)  +  ^2(0$)  A-j  ±  2  Ai  Ag  sin(^i  -  ij^)  +  Ag 


1/2 


cos(</?  +  a  +  i|>+) 


(41) 


where 


P  ~  Aj/A^  » 


a  s  tan  [Gg(8 ^)/Gi  (9 ^)] 


ip+  -  tan 


,1/A1  sin  ^  1  ±  A2  cos^ 
\A1  cos  ^  ^  +  sin  ^ 


(42a) 

(42b) 

(42c) 
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Unfortunately,  there  are  too  many  parameters  involved  in  these  equations  to 
obtain  any  physical  'nsight  into  the  perform? nee  of  the  detector  in  the  general 
case.  We  have  had  to  resort  to  using  a  simulation  to  evaluate  this  general  case. 
We  can  ottain  some  useful  analytical  results  by  limiting  our  interest  to  real 
antenna  patterns.  This  will  be  valid  for  ma inbeam  interference  and  also  for 
interference  located  at  the  near-in  sidelobes.  Therefore,  we  assume  that  G. (sT) 
are  real.  In  this  case,  (41)  becomes 


where 


K<es>  +  Gl(es> 

+  p2  G 

+  2p 

,  ?  ’/2 
G,(8S)  +g|(6s) 

Gl(6$) 

G-j  (flj) 

tan"^ 

w 

:  +  giK> 

rw 

G^e'j) 

Ui 


1/2 


(44) 


It  will  be  useful  to  use  the  notion  of  the  monopulse  function  which  is  defined 
as 


E(e) 


G^e) 

G-j  ( 9 ) 


(45) 


In  most  practical  monopulse  systems,  the  monopulse  function  is  linear  over  the 
extent  of  the  main  beam.  Since  the  target  of  interest  will  always  be  within 
the  mainbeam  then 
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(46) 


E(6S)  = 


where  eg  is  the  3  dB  beamwidth  and  k  is  a  standard  parameter  that  arises  in  the 
analyses  of  monopulse  systems.  Typicall>,  1  =  k  =  2  with  k  -  1.5  being  a  reason¬ 
able  value  [8].  -  Whereas,  the  target  can  always  be  assumed  to  lie  within  the 
mainbeam,  the  interference  signal  can  originate  from  the  mainbeam  or  from  a  side- 
lobe.  Howe/er,  it  Is  convenient  to  define  an  equivalent  interference  azimuth, 

6j  as 

~*E(aT)/k  .  (47) 

0B  1 


We  note  that  e  =  0j  when  is  within  the  mainbeam  of  the  antenna.  In  Fig.  2a 
we  have  plotted  9j  vs  6j  using  a  typical  monopulse  function  derived  for  a  4° 
beamwidth  antenna  with  -20  dB  sidelobes,  where  errors  in  the  amplitude  and  phasp 
tapers  render  the  antenna  pattern  complex.  Therefore,  the  results  can  be  expected 
to  give  some  indication  of  performance  even  in  the  more  general  case.  Assuming 
Gj  is  uniformly  distributed  in  ( -tt ,tt) ,  Fig.  2b  gives  the  probability  distribution 
function  of  0j.  This  shows  that  0j  will  lie  within  ±  2  beamwidths  most  of  the 
time  and  hence  the  cases  of  sidelobe  and  mainbeam  interference  can  be  treated 
simultaneously. 

Using  (46)  we  can  write  E(6j)  =  k  Bj/Bg,  and  then  absorb  eg  into  our 
definition  of  e  so  that  all  of  the  azimuth  variables  can  be  expressed  in  3  dB 
beamwidths.  Using  these  relations  (42)  and  (43)  can  be  written  as 


d*  =  4  Gi 


<8S>  (’  +  k2  9s) 


1  +  p 


2/'  *  k  4 


[1  +  k  Q 


ry/+  \ 


(1  +  k2  e2^2 

— FT 


COS(c,o  +  B) 

(48) 
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mmmm 


where 


g  =  tan 


-1 


k(es  -  0j ) 

l  +  k2  es  9j 


Aj  ( 0  •£ ) 

po  =  As  &1  res ) 


(49) 


(50) 


The  quantity  Pq  represents  the  interference-to-signal  ratio  as  measured  at  the 
output  of  the  antenna  terminals.  We  note  that  by  using  (47)  we  are  able  to 
derive  an  expression  which  describes  the  case  of  mainbeam  and  sidelobe  inter¬ 
ference  simultaneously. 

a)  Basic  Properties  of  the  Detector 

Using  the  above  expression  we  can  obtain  some  interesting  properties 
of  the  interference  detector.  First  we  note  that  if  the  relative  phase  between 
the  two  signals,  (/>,  is  0  or  w,  then  d+  =  d_  and  from  (39)  we  see  that  the  detec¬ 
tion  probability  reduces  to  the  false  alarm  probability.  This  is  an  unfortunate 
property  since  McAulay  [1]  has  shown  that  interference  causes  the  worst  azimuth 
errors  at  the  out-of-phase  condition.  Secondly,  we  note  that  if  8<.  =  0j  =  0j  then 
S  =  0  and  d+  =  d_.  Again,  the  detection  probability  is  negligible.  This  is  a 
reasonable  behavior  for  the  detector  to  exhibit  since  if  targets  are  at  the 
same  azimuth,  they  will  not  cause  any  azimuthal  error  except  that  due  to  fading  [1]. 

Another  interesting  analytical  result  can  be  obtained  in  those  cases  where 
the  interference  has  been  completely  overpowered  by  the  target.  In  this  case, 
pQ  «  1  and  we  can  neglect  the  second  term  in  (48)  and  use  the  Binomial  Expansion 
to  take  the  square  root.  This  gives 
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Using  (49)  it  can  then  be  shown  that 


d 


+ 


”  Aj  G-j  ( 0j  ; 

2a  0 


k(  65  -  Qj) 

0  +  k2  e^)1/2 


sin  <p 


(52) 


This  result  shows  that  when  the  target  completely  overpowers  the  interferes  the 
ability  of  the  receiver  to  detect  interference  depends  on  the  interference-to- 
noise  ratio  (INR).  At  first  glance  this  is  a  somewhat  puzzling  result  since 
McAulay  [1]  has  shown  that  the  azimuth  accuracy  depends  on  the  signal-to-inter- 
ference  ratio  (SIR)  such  that  if  the  SIR  is  large,  the  azimuth  estimate  is  un¬ 
affected  by  interference.  Equation  (52)  indicates  that  if  the  INR  is  also  large 
then  the  interference  detector  would  ring.  From  a  data  editing  point  of  view 
this  would  not  be  a  desirable  property.  However,  the  detector  is  not  really 
testing  for  the  presence  of  interference  since  there  is  no  inherent  distinction 
between  interference  and  target.  Rather  it  is  testing  for  the  presence  of  more 
than  one  signal.  If  the  SIR  is  very  small,  then  (48)  reduces  to 


d+  -  d  Ar  G-.  (9$)  k(eT  -  9c) 

—5 — 1  - sin  V 

2a  0  (1  +  kZ  9Z)1/Z 


(53) 


In  this  case,  a  large  azimuth  error  would  result  since  the  azimuth  of  the  strong 
interferer  would  be  estimated.  However,  this  situation  would  be  flagged  by  the 
detector  provided  the  target  SNR  were  large  enough.  Therefore  the  detector 


performance  in  these  extreme  cases  is  intuitively  satisfying,  although  the 
results  indicate  that  some  provision  may  have  to  be  made  for  reducing  the  de¬ 
tections  due  to  low  level  interference.  We  will  discuss  this  point  further  in 
a  later  section. 

To  obtain  a  more  complete  understanding  of  the  detector  performance  we 

need  to  evaluate  (48)  and  (49)  and  use  these  with  (39)  and  (37)  to  obtain  the 

detection  and  false  alarm  probabilities.  The  results  we  obtain  are  based  on  the 

Gaussian  approximation  to  the  Rician  densities.  From  (28)  we  see  that  this  will 

2  2  > 

be  a  reasonable  assumption  provided  d+/2a  =  10.  From  (48)  we  see  that  the 
2 

smallest  value  of  d+  occurs  when  ip  +  0  =  tt.  Then 


2  a 


7- 


AS  Gl(eS> 


2  „2> 


2  a 


'  (1  +  k  0g) 


h  ♦  k2  e2\'/2 

- FT 

,1  ♦  k2  e2y 


=  10 


(54) 


By  requiring  that  the  second  inequality  hold,  we  have  a  conservative  but  sufficient 
condition  to  guarantee  that  (28)  will  hold.  Therefore,  for  a  post-detection  SNR 
greater  than  20  dB,  t.nis  inequality  will  be  satisfied  for  all  signal-to-inter- 
ference  ratios  except  those  in  the  region  from  -3  dB  to  +  3  dB.  Since  this  is 
a  conservative  assumption  and  since  we  would  not  expect  the  detection  probability 
to  depend  on  the  tails  of  the  pdf,  it  is  reasonable  to  expect  that  the  Gaussian 
approximation  will  adequately  describe  the  pfcfformance  over  an  even  smaller 
range  of  SIR  values. 

Over  the  duration  of  any  one  DABS  reply,  the  parameters  A<.,  Aj ,  6$,  6j  are 
fixed.  For  mainbeam  multipath,  (p  will  also  be  constant  from  chip  to  chip.  For 
the  case  of  ATCRBS  interference  <pwill  change  randomly  from  chip  to  chip,  and 
in  fact  may  also  change  between  samples  within  a  chip  due  to  the  frequency  offset 


that  can  be  expected  between  the  OABS  and  ATCRBS  transponders.  In  addition,  we 
have  found  that  the  dependence  on  <p  is  crucial  to  the  understanding  of  the  data 
editing  concept.  Therefore,  we  have  chosen  to  evaluate  the  detection  probability 
as  a  function  of  <p. 

b)  Performance  Based  on  a  Single  Sample 

We  have  evaluated  the  performance  of  the  receiver  for  some  cases 
of  interest.  Our  basic  parameter  values  were  taken  to  be  the  following: 

6S  =  0  (target  on  bores ight)  , 

Aj  Gj(0) 

SNR  =  -  y~~  =  20  dB  , 

2  a 

=  0.5  (interference  at  3  dB  point)  ,  (55) 

k  =  1.5  , 

G^e)  =  l  -  1.17  e2  . 

All  we  need  do  now  is  specify  a  false  alarm  probability,  compute  XpA  from  (37), 
compute  <\J 2c  from  (48),  (29)  and  use  this  to  compute  the  detection  probability 
in  (39).  In  Figure  3  we  have  plotted  the  normalized  mean  value  of  the  detection 
statistic  (d+  -  d_)/2a  for  several  values  of  the  ISR. 

Since  in  the  DABS  direction  finding  (DF)  problem  there  will  be  many  bits 
available  for  generating  the  azimuth  estimate,  a  higher  false  alarm  rate  can  be 
tolerated  in  order  to  correctly  detect  interference  samples.  Therefore,  we 
allowed  a  false  alarm  rate  of  2  samples  in  100,  PFA  =  0  .02,  and  then  computed  the 
detection  probability  fo**  various  values  of  the  ISR. 
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In  Figure  4  we  have  plotted  the  miss  probability  P^(c>)  =  1  -  P^Gp)  as  a  function 
of  the  relative  phase.  As  we  expected,  the  detector  misses  interference  with 
probability  one  when  the  in-phase  and  out-of-phase  conditions  exist.  In  the  case 
of  ATCRBS  interference  the  relative  phase  is  independent  and  uniformly  distributed 
from  bit  to  bit.  Then  a  measure  of  performance  of  the  detector  is  the  average 
miss  probability 


r  2* 

PM  =  h  pni<p)  d(p 

Jo 


(56) 


In  Figure  4  we  have  indicated  the  average  miss  probability  as  a  function  of  ISR 
for  the  0.02  false  alarm  probability  case.  j 

c)  The  Effect  of  Frequency  Offset  ! 

The  preceding  results  give  the  detection  performance  when  a 
decision  must  be  made  on  the  basis  of  a  single  sample.  In  the  DABS  context, 
there  will  be  several  samples  available  per  chip  for  interference  detection  and 
direction  finding.  Unfortunately,  if  the  interference  is  multipath,  the  phase 

j 

will  not  change  significantly  from  sample  to  sample,  or  even  from  reply  to  reply,  I 

hence  the  detector's  performance  will  be  essentially  the  same  as  that  described  | 

in  the  last  section.  Hence  there  may  be  situations  when  azimuthal  multipath  | 

will  be  present,  but  will  not  be  detected  by  the  interference  flag.  | 

| 

In  the  case  of  ATCRBS  interference,  however,  there  will  be  a  frequency  ,  j 

i 

offset  between  the  DABS  and  ATCRBS  transponders  that  can  cause  the  instantaneous  ] 

3 

4 

phase  to  change  from  sample  to  sample.  For  example,  if  &f  denotes  the  frequency  j 

\ 

difference  between  the  two  transponders  and  if  fg  represents  the  rate  at  which  j 

the  DABS  waveform  is  to  be  sampled,  then  from  sample  to  sample  the  phase  increases  j 

i 

} 
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RELATIVE  PHASE  BETWEEN  TARGET  AND  INTERFERENCE  SIGNALS  (rad) 

Fig,  4.  Probability  interference  is  not  detected. 
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by  2ir  Af/f$.  Since  the  sampling  is  to  be  done  at  a  10  MHz  rate,  then  a  0.5  MHz 
frequency  offset*  would  lead  to  a  0.31  radian  phase  shift.  From  Figure  4  we  see 
that  the  detector  misses  the  interference  when  the  phase  is  within  a  1  radian 
interval  about  0  or  n.  Therefore,  with  two  or  three  additional  samples  we  can 
expect  the  detector's  performance  to  improve  significantly.  We  can  make  these 
statements  quantitative  by  considering  the  detection  of  interference  using  N 
samples  per  chip.  Our  strategy  is  to  declare  interference  present  if  the  detection 
threshold  is  crossed  for  at  least  one  of  the  N  samples.  Then  to  yield  a  miss,  the 


detector  must  fail  to  detect  on  every  sample.  If  cp  denotes  the  relative  phase 
at  the  time  of  the  first  sample,  then  at  the  nth  sample  it  is  <p  +  (n-1 )  Ac?.  The 


single  sample  miss  probability  at  this  phase  is  denoted  PM  \jp  +  (n-1)  Ac?].  Then 


the  miss  probability  after  N  samples  is 


N 

PM  M  =rTPM  (n-1)  A^]  .  (57) 

N  W  1 

The  product  rule  applies  because  the  noise  samples  are  independent  from  sample 
to  sample.  If  we  let  T  denote  the  chip  width  and  f  the  sampling  rate,  then  the 
number  of  samples  per  chip  is  N  =  f  T,  and  the  phase  shift  is  Ac?  =  2ir  Af/f$.  For 
the  DABS  application  we  could  expect  T  =  0.483  usee  and  f  =  10  MHz.  Therefore 
3  samples  per  chip  represents  a  conservative  evaluation.  Using  this  value  for 
N,  we  plot  PM  as  a  function  of  frequency  offset  for  several  values  of  the  ISR 
and  an  0.02  false  alarm  rate.  This  is  shown  in  Figure  5  and  demonstrates  the 
significant  improvement  in  detection  performance  that  can  be  expected  as  a  result 


distribution  of  frequency  offset  has  been  measured  by  G.  Colby  and  E.  Crocker 
and  is  documented  in  reference  [9]. 
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AVERAGE  DETECTION  PROBABILITY 


0  0.5  1  1.5 

FREQUENCY  OFFSET  (MHz) 


Fig.  5.  Effect  of  frequency  offset 


of  the  instantaneous  phase  change  from  sample  to  sample  that  results  from  the 
frequency  offsets  between  a  DABS  and  an  ATCRBS  transponder. 
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VI.  CONCLUSIONS 

The  analysis  of  the  performance  of  the  maximum  likelihood  interference 
detector  was  greatly  simplified  by  recognizing  that  the  interference  statistic 
could  be  well -approximated  by  a  Gaussian  random  variable.  A  conservative  condition 
on  when  the  approximation  was  valid  was  found.  For  example,  at  a  20  dB  SNR,  the 
SIR  would  have  to  be  within  +  2  dB  for  the  approximation  not  to  be  valid. 

False  alarm  and  detection  probabilities  were  calculated  in  detail  for  the 
cases  in  which  the  antenna  patterns  were  real.  Although  formulae  for  the  complex 
antenna  pattern  case  were  derived,  it  was  not  possible  to  obtain  simple  analy¬ 
tically  useful  results  from  them.  We  then  restricted  our  attention  to  the  case 
of  real  antenna  patterns  and  obtained  expressions  from  which  it  was  possible  to 
draw  some  useful  conclusions.  It  was  found  that  when  the  relative  phase  was  C 
or  tt,  the  receiver  would  fail  to  detect  the  interference  with  probability  one. 

It  was  also  noted  that  the  detectability  improved  as  the  azimuthal  separation  of 
the  two  signal  sources  increased.  Furthermore,  it  was  found  that  detectability 
depended  on  the  signal-to-noise  ratio  of  the  weaker  of  the  two  signals.  If  this 
SNR  was  large  enough,  good  detection  was  obtained  no  matter  how  large  the  other 
signal  became. 

Specific  results  for  a  mainbeam  interferer  were  given  and  it  was  found  that 
if  only  a  single  sample  were  used  for  detection  that  the  average  detection  prob¬ 
ability  was  approximately  0.8  for  an  0.02  false  alarm  rate.  This  poor  perform¬ 
ance  was  due  to  the  fact  that  misses  were  guaranteed  when  the  phase  differences 
were  0  or  tt. 
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In  practice  there  will  be  several  samples  per  chip  available  for  inter¬ 
ference  detection  which  can  result  in  meaningful  improvements  in  performance 
provided  there  is  an  instantaneous  phase  change  from  sample  to  sample.  Un¬ 
fortunately,  if  the  interference  is  due  to  multipath,  there  will  be  no  signifi¬ 
cant  change  in  phase  for  several  seconds  duration,  hence  multiple  samples  cannot 
be  expected  tc  improve  the  performance  in  this  case.  When  the  interference  is 
due  to  ATCRBS,  however,  the  probable  frequency  offset  between  the  DABS  and  ATCRBS 
transponders  will  cause  the  phase  to  change  from  sample  to  sample.  We  have  ex¬ 
amined  this  case  in  detail  and  found  that  significant  improvements  in  target 
detection  performance  can  be  expected. 
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APPENDIX 

DERIVATION  OF  THE  APPROXIMATE  PDF 


Let  us  focus  our  attention  on  the  Rician  pdf 


P(“)  *  T  exp(-  4^)  ’ol^) 


The  well-known  asymptotic  approximation  for  the  I0(*)  Bessel  Function  is 


(A-l) 


’o^)'  MH)  7  cxp(!7) 


(A -2) 


This  is  quite  an  accurate  approximation  for  ud /<j  >_  3.  Then  the  pdf  becomes 


PW*  <$yt -T=«>  ‘  1*7$ 

V2"0 


(A-3) 


We  are  most  interested  in  the  values  of  u  about  d.  Let  us  write 


u  =  d  +  6  u 


(A-4) 


<1  +  r) 


-t— -  exp - 

\i%  n  2  '  20 

V2"  ° 


(-$ 


=  (1  +  \  4r-  +  . . . )  -- ~r-  exp 


(A-5) 


Therefore,  u  is  well  approximated  by  a  Gaussian  random  variable  with  mean  d  and 
variance  2  a2  provided  ud/o2  =  3  and  ^  «  1.  When  the  Gaussian  approxi¬ 
mation  is  valid,  95%  of  the  pdf  lies  between  the  2-sigma  limits.  Therefore,  if 
we  require  that 


1.  2o<  1 

2  d 


(A-6) 


then  j  «  1  "most  of  the  time."  As  a  practical  matter  we  take 


o  <  1 


/io 


(A-7) 


as  the  criterion  for  which  we  can  neglect  first  and  higher  order  terms  of  Su 
that  appear  in  (A-5).  We  must  yet  determine  whether  or  not  (A-7)  is  sufficient 
to  validate  the  IQ(*)  approximation  that  led  to  (A-5)  in  the  first  place. 
However,  u  =  d  -  2a  at  least  97.5%  of  the  time,  hence 


H_|+d^  0  .  20)  >  ^  .  JL_)  =  2.7  . 

~  -  a  '  /to 7 


(A-8) 


Therefore,  the  condition 


d2> 


10 


(A.-9) 


is  sufficient  to  guarantee  that  the  Rician  pdf  for  u  will  be  well  approximated 
by  the  Gaussian  density  of  mean  d  and  variance  rr‘\ 
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